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ABSTRACT 

The  question  of  optimizing  nozzle  contours  for  micro-thrust 
rockets  led  to  the  design,    construction,    and  testing  of  a  low-density 
gas  dynamic  facility.      The  primary  objective  was  to  investigate  the 
mass  flow  rates  of  a  gas  through  various  profiles  in  the  slip  and 
transition  flow  regimes  at  high  pressure  ratios. 

An  initial  test  was  conducted  with  an  orifice  as  the  test  pro- 
file.     The  results  showed  that  the  facility  can  be  used  to  investigate 
mass  flow  rates  from  the  threshold  of  the  free-molecule,    through 
the  transition  and  slip,    to  the  continuum  regimes.      These  results 
compare  favorably  with  those  of  two  previous  investigators,    and 
asymptotically  approach  the  theoretical  continuum  and  free-molecule 
limits.     The  ratio  of  mass  flow  rate  to  theoretical  free-molecule 
mass  flow  rate  is  shown  to  transition  smoothly  from  one  theoretical 
limit  to  the  other.     A  local  maximum  may  occur  in  this  ratio  in  the 
slip  regime,  and  the  attainment  of  the  theoretical  free-molecule  limit 
appears  to  occur  more  slowly  than  expected. 
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a  molecular  collision  diameter 

c  mean  molecular  speed 

d  aperture  diameter 

f  transfer  ratio  of  momentum 

g  gravitational  constant 

h  height  difference 

k  constant  '  <  ' 
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p  pressure 
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F  collision  frequency  per  mole 

H  perimeter 
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M  molecular  weight 

Ma  Mach  number 

N  number  of  molecules  striking  a  unit  area  of  wall  per  unit 
time 

P  pressure 
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Q  volumetric  flow  rate 

R  specific  gas  constant 

Re  Reynolds  number 

T  temperature 

V  volume 

W  characteristic  velocity 
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Q  coefficient  of  slip 

0  spherical  elevation  angle  or  conical  half  angle 
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I.     INTRODUCTION 
The  flow  of  gases  can  be  roughly  divided  into  four  regimes 
characterized  by  the  Knudsen  number  (Kn),    the  ratio  of  molecular 
mean  free  path  to  a  characteristic  length  in  the  flow  field.      These 
four  regimes  are: 

continuum  flow  ~  Kn  <  0.  01 

slip  flow  ~    0.  01  <  Kn  <  0.  1 

transition  flow  ~    0.  1     <  Kn  <  1 

free-molecule  flow  ~    5  <  Kn 

The  Navier-Stokes  equations  with  appropriate  boundary 
conditions  describe  the  flows  in  the  continuum  and  slip  flow  regimes 
while  the  kinetic  theory  of  gases  must  be  used  in  the  transition  and 
free-molecular  regions. 

For  flows  through  orifices,   nozzles, and  short  tubes  well 
into  the  continuum  regime,    the  boundary  layers  on  the  walls  usually 
represent  only  a  very  small  portion  of  the  flow  field.     These  boundary 
layers  become  more  and  more  dominant  as  the  Knudsen  number  in-    - 
creases,    until  they  occupy  the  entire  flow  field  at  large  Knudsen 
numbers.     In  the  slip  regime,    the  boundary  layers  no  longer  have 
zero  velocity  at  the  walls.     In  free-molecule  flow, collisions  between 
gas  molecules  can  be  neglected  compared  with  collisions  between  gas 
particles  and  the  body.     In  such  cases,    the  fluxes  of  incident  and  re- 
flected particles  can  be  treated  separately.     In  the  transition  regime 
where  the  mean  free  path  and  the  characteristic  dimension  of  the 
body  are  of  the  same  order  of  magnitude,   both  intermolecular  colli- 
sions and  molecule -surface  interactions  must  be  considered.     In 
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order  to  obtain  experimental  data  for  a  better  understanding  of  the 
transition  from  continuum  to  free-molecule  flow,  the  present  low- 
density  gas  dynamic  facility  was  designed  and  built. 

The  initial  interest  in  this  area  of  low-density  gas  flow 
was  generated  by  the  problem  of  nozzle  contours  for  micro-thrust 
rocket  motors.     In  these  rockets  with  extremely  low  thrust,    the  ex- 
haust gases  are  in  the  flow  regime  between  continuum  and  free- 
molecule.      The  optimization  of  these  nozzle  contours  depends  upon 
having  some  basic  knowledge  of  gas  flow  in  the  slip  and  transition 
regimes.     It  is  hoped  that  this  low-density  gas  dynamic  facility  will 
supply  some  basic  experimental  data  about  slip  and  transition  flow 
through  various  contours. 
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II.     APPARATUS 
General  Description 

The  low-density  gas  dynamic  facility  is  a  continuous  flow, 
open-circuit  type  designed  to  operate  with  stagnation  chamber  pres- 
sures of  from    1  [l  to  \  atm   and  exhaust  chamber  pressures  to  give 
at  least  a  100:1  pressure  ratio  across  the  test  section. 

Figure  1  is  a  photograph  of  the  facility;  subsystems  are 
shown  schematically  in  Figures  2  -  4  .     The  pumping  system  consists 
of  a  300-cfm  rotary  pump,    Stokes  Model  41 2- H,    in  parallel  with  a 
2-in.    metal-oil  fractionating  diffusion  pump,    CVC  Model  PMCS-2C, 
backed  by  a  7-cfm  rotary  pump,    Cenco  Hypervac  25  (Figure  4).     The 
pumping  characteristics  as  supplied  by  the  manufacturers  are  shown 
in  Figures  5  and  6. 

The  vacuum  vessel  (Figure  3)  is  a  cylindrical  steel  cham- 
ber of  about  40-cm  internal  diameter  and  90-cm  length.     The  tank  is 
divided  into  the  stagnation  and  exhaust  chambers  by  a  central  bulk- 
head which  contains  a  cutout  for  mounting  a  machined  test  profile 
plate.     For  rapid  pump-down,    a  bypass  is  incorporated  around  the 
test  section.     The  tank  has  provisions  for  three  stagnation- chamber 
pressure  taps.     The  exhaust  chamber  pressure  line  enters  the 
chamber  through  the  4-in.    outlet  and  extends  about  half  way  into  the 
chamber.      The  tip  of  the  pressure  line  is  off  to     the  side  of,    and 
turned  at  a  right  angle  to,    a  direct  line  between  the  test  section  and 
•  the  outlet.     Access  to  the  test  section  is  gained  through  an  O- ring- 
sealed  plate  in  the  inlet  side  of  the  stagnation  chamber.     This  plate 
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incorporates  the  inlet  line  and  one  pressure  tap. 

Test  Gas  Feed  and  Measurement 

The  test  gas  is  fed  from  a  standard  gas  bottle  through 
first,    a  25-  to  650-psig  high-pressure  regulator,    Matheson  Model  2, 
and  then,    a  3-  to  15-in.    water-column  low  pressure  regulator, 
Matheson  Model  70  B,    before  passing  to  a  liquid-nitrogen  cold  trap. 
From  the  cold  trap,    the  gas  flows  through  a  constant-temperature 
water  bath  to  the  flow  measurement  system.     At  the  inlet  to  the  flow 
measurement  system,    there  are  both  pressure  and  temperature  taps. 
The  pressure  line  is  connected  to  a  0-  to  15-psia  dial  gauge,    Burton 
Model  2311-001  and  a  0-  to  800-mm  mercury  manometer,    Wallace 
and  Tiernan  Model  FA  135,    and  the  temperature  tap  houses  a  pre- 
cision mercury  thermometer,   ASTM-Saybolt  Viscosity  Thermometer 
17F  (66     to  80     F).     The  pressure  gauge  can  be  read  to  within  ±  0.  01 
psia  and  is  accurate  to  within  \  per  cent  full  scale.     The  thermome- 
ter can  be  read  to  within  ±  0.  05     F. 

The  flow  measurement  system  consists  of  two  subsystems 
of  different  flow  measuring  devices  (Figure  2).     The  first  is  a  paral- 
lel array  of  four  rotameters,    Brooks  Models  R-2-15-AA  with  stain- 
less steel  float,    R-2-25-D,    R-2-25-A,    and  R-2-25-B  with  sapphire 
floats.     These  four  rotameters  are  capable  of  measuring  helium 
flow  over  the  ranges  0.  4  to  4.  3,    1.  0  to  10.  9,    4.  0  to  36,    and  10  to 
108  cc/sec  at  760  mm  Hg  and  70     F.     The  rotameters  were  calibrated 
by  the  author  after  they  were  installed  in  the  facility  by  the  water 
displacement  method  at  atmospheric  conditions.     Each  calibration 
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curve  is  based  on  at  least  20  data  points  that  are  reproducible  to 
within  ±  1  per  cent  from  30  to  100  per  cent  of  full-scale  reading. 
The  overlapping  range  of  adjacent  rotameters  permits  each  meter  to 
be  used  in  its  upper,    more  accurate,    ran^e  (Figure  7). 
The  governing  equation  for  a  rotameter  is: 

2gVf(pf-pw) 


where 


Q    =    A    C 
w 


2 


A,p 


fw 


(1) 


Q    =    volumetric  rate  of  flow, 
V-    =    volume  of  float, 
g    =    gravitational  constant, 
p,    =    float  density, 

p       =    fluid  density 

rw  ' 

A-    =    area  of  float, 

C    =    discharge  coefficient, 

_.  2      2 

A       =    area  of  annular  orifice,     =    -r[(D+by)    -d   ], 

D    =    effective  diameter  of  tube  depending  on  position  of 
float, 

b    =    change  in  tube  diameter  per  unit  change  in  height, 

d    =    maximum  diameter  of  float,    and 

y    =    height  of  float  above  zero  position. 
The  value  of   C   is  a  weak  function  of  Reynolds  number,    but  through 
careful  design,    the  effect  is  largely  eliminated  for  certain  ranges  of 
Reynolds  number.     Thus,    when  the  rotameter  is  calibrated  and  used 
for  a  single  fluid  at  approximately  the  same  conditions  of  temperature 
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and  pressure,     C   may  be  taken  as  a  constant. 

For  a  gas    p      is  negligible  compared  with   p.  ,    so  equation 
(1)  becomes 

2gVfpf 


Q    =    A    C 

w 


ii 


A,p 
frw 


(2) 


Therefore,   for  any  given  position  of  the  float    Q  =  kp     2  where   k  is 
a  constant   . 

Once  the  meter  is  calibrated  at  known  conditions  of  tem- 
perature and  pressure  for  a  given  gas,    the  reading  can  be  interpreted 
at  other  (not  greatly  different)  conditions  for  the  same  gas  by  know- 
ing the  gas  densities  or  equation  of  state: 

rw  c 

where  the  subscript   c    stands  for  the  calibration  condition.     The 
mass  flow  rate  (m)  therefore  can  easily  be  determined: 

m    =    Qp       =    Q   (p       p    }z    .  (4) 

rw  c  rwcrw 

The  second  flow  measurement  subsystem  is  a  Brooks 
Vol-U-Meter.     This  is  a  precision-bore,   vertical  borosilicate  glass 
tube  containing  a  polyvinyl  chloride  floating  piston  having  a  mercury 
O-ring  seal.     The  usable  length  of  the  tube  is  30  in.    and  the  internal 
volume  is  indicated  every  1  cc  from  0  to  25  cc.      The  precision  of 
the  instrument  is  0.  2  per  cent  of  indicated  volume.     3y  measuring 
the  time  required  for  the  piston  to  rise  through  a  given  displace- 
ment and  knowing  the  output  pressure  and  temperature,    the  mass 
flow  rate  can  be  calculated.     Volume  flow  rates  from  about  0.  003 


-7- 
to  2.  5  cc/sec  can  be  measured,    the  lower  end  being  limited  only  by 
the  patience  of  the  operator  and  the  overall  stability  of  the  system, 
and  the  upper  end  by  the  reaction  time  of  the  operator  and  the  diffi- 
culty in  gauging  when  the  piston  passes  a  given  volume  mark. 

It  was  found  that  the  test  gas  stays  at  room  temperature 
throughout  the  facility  after  it  has  passed  through  the  water  bath  at 
room  temperature.      Thus,    the  temperature  at  the  input  to  the  flow 
measuring  system  is  taken  as  the  temperature  in  the  Vol-U-Meter. 
A  pressure  tap  is  installed  in  the  output  line  of  the  Vol-U-Meter. 
This  pressure  is  indicated  on  a  0  to  25-  and  25  to  50-psig  dual- 
range  dial  gauge,    Wallace  and  Tie  man  Model  FA2331H,    and  is 
taken  as  the  pressure  in  the  Vol-U-Meter.     It  was  found  that  the 
pressure  drop  across  the  mercury-sealed  piston  was  about  0.  08  psi, 
or  less  than  1  per  cent  of  the  input  pressure. 

During  runs  in  which  the  Vol-U-Meter  is  to  be  used,    the 
Vol-U-Meter  is  bypassed  until  a  steady-state  condition  is  attained. 
Once  steady-state  is  established,   the  Vol-U-Meter  bypass  is  closed 
and  the  test  gas  is  drawn  from  the  Vol-U-Meter  until  the  flow  rate 
has  been  determined.     Once  this  has  been  accomplished,   the  bypass 
is  again  opened.     The  steady-state  conditions  are  checked  before, 
during,    and  after  the  volumetric  measurement  in  order  to  be  certain 
that  the  stagnation  chamber  remains  at  the  same  steady-state  condi- 
tion. 

From  the  flow  measurement  system,    the  test  gas  is  fed  to 
the  stagnation  chamber  through  a  variable   leak,    Granville -Phillips 
Series  203.     This  variable  leak  is  used  to  control  the  test-gas  flow 
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rate,    and  thus  the  pressure  in  the  stagnation  chamber.      The  con- 
ductance of  the  variable  leak  is  continuously  variable  from  approxi- 
mately 100  to  10  cc/sec  with  atmospheric  pressure  on  the  inlet 

side  and  vacuum  on  the  outlet.      Fully-closed  conductance  is  stated 
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by  the  manufacturer  to  be  less  than  10  cc/sec.     Twenty-seven 

full  turns  of  the  driver  handle  provides  precise  control  over  the 
entire  conductance  range.     A  counter,    graduated  in  units  correspond- 
ing to  1/10  turn  of  the  handle,    is  helpful  in  setting  approximate  con- 
ductances.    At  the  higher  flow  rates  and  stagnation  chamber  pres- 
sures,   a  given  setting  gives  reproducible  performance.     At  the  lower 
settings,   hysteresis  losses  found  in  all  devices  with  metal-to-metal 
seals  cause  slight  variations  in  performance  at  any  given  setting. 
Figure  8  shows  approximate  counter  setting  versus  stagnation  cham- 
ber pressure  over  the  range  of  0.  01  to  200  mm  Hg. 

Test  Profile  Plate 

The  test  profiles  are  machined  in  steel  plates  which  are 
fitted  with  O-ring  seals  into  a  cut-out  in  the  dividing  bulkhead  of  the 
vacuum  vessel.     The  stagnation  chamber  side  of  the  plates  is  flush 
with  the  surface  of  the  dividing  bulkhead  so  that  no  undue  entrance 
effects  are  present.     The  minimum  diameter  of  the  test  profile  is 
approximately  1  mm.      Thus,    the  ratio  of  tank  diameter  to  profile 
diameter  is  about  400  to  1.     For  all  practical  purposes,    the  stag- 
nation and  exhaust  chambers  are  of  infinite  size  in  comparison  with 
the  test  profile.     Therefore,    the  effects  of  the  tank  wa^is  on  the  flow 
field  through  the  test  profile  are  negligible. 
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The  first  profile  tested  was  a  sharp-edged  orifice   with  a 
diameter  of  0.  039  in.    and  an  edge  thickness  of  0.  0035  in.    (Figure  9). 
Subsequent  profiles  will  be  short  tubes  and  conical  profiles  of  vari- 
ous entrance  and  exit  half  angles. 

Pump  System 

The  4-in.    outlet  of  the  exhaust  chamber  is  connected  to  a 
welded  aluminum  tee.      The  tee  has  4-  and  2  l/4-in.    flanged  outlets. 
Both  the  outlets  are  fitted  with  handwheel-ope rated  gate  valves, 
Temescal  Series  5000.     From  the  4-in.   valve  the  flow  drops   12  in. 
in  a  flanged  Zee  connector  which  serves  as  a  baffle  to  prevent  back- 
flow  of  oil  from  the  Stokes  pump.     The  Stokes  pump  is  vibration  iso- 
lated  from  the  Zee  connector  by  means  of  a  6-in.    stainless  steel  bel- 
lows.    The  exhaust  of  the  Stokes  pump  is  vented  outside  the  building 
(Figure  1). 

The  2  1 /4-in. Temescal  valve  is  connected  by  O-ring  seals 
to  the  input  flange  of  the  CVC  diffusion  pump.      The  exhaust  of  this 
pump  is  fed  to  the  Cenco  rotary  backing  pump  by  means  of  1-in. 
diameter  copper  tubing  that  is  isolated  from  the  diffusion  pump  by 
means  of  a  short  piece  of  Tygon  tubing  (Figure  1). 

This  pumping  arrangement  permits  rapid  pump-down  of  ' 
the  system  to  the  micron  range  by  use  of  the  Stokes  pump.     The 
pump-down  time  to  10  |j   is  approximately  one  minute  (Figure  10). 
The  diffusion  pump  can  be  isolated  from  the  system  and  warmed  to 
operating  conditions  while  the  Stokes  pump  is  exhausting  the  system. 
The  switch-over  to  the  diffusion  pump  is  accomplished  without  any 
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increase  in  system  pressure  by  simply  opening  the  2  1/4-in.    Tem- 
escal  valve  and  closing  the  4-in.   valve.     It  has  been  :he  practice  to 
keep  the  diffusion  pump  in  continuous  operation  and  simply  isolating 
it  from  the  system  when  the  flow  rate  and  exhaust  chamber  pressures 
are  such  as  to  dictate  the  use  of  the  Stokes  pump  as  the  primary 
vacuum  source. 

Vacuum  Pressure  Measurement 

The  primary  means  of  measuring  vacuum  pressure  is  a 
bank  of  two  McL.eod  gauges,    made  by  Greiher  Glassblowing  Labora- 
tory,   a  0-  to  800-mm  mercury  manometer,    Wallace  and  Tiernan 
Model  FA  135,    and  a  0-  to  50-mm  dial  gauge,    Wallace  and  Tiernan 
Model  FA  160.     One  McLeod  gauge  is  a  non-linear  gauge  that  meas- 
ures the  range  0  to  125  |j   on  a  230-mm  scale  and  that  has  a  capture 
volume  of  330.  6  cc.     The  second,    linear  gauge,   has  three  500-mm 
scales  that  measure  the  ranges  0  to  1,    0  to  10,    and  0  to  100  mm 
with  a  capture  volume  of  250.  5  cc. 

In  the  non-linear  gauge,    a  known  volume  of  sample  gas  at 
the  unknown  pressure  is  captured  and  then  compressed  in  a  1-mm 
capillary  that  is  sealed  at  its  upper  end  until  a  reference  column  of 
mercury  in  a  similar  capillary  is  even  with  the  end  of  the  sealed 
capillary.     The  difference  in  the  heights  of  the  mercury  columns  in 
the  two  capillaries  is  then  a  measure  of  the  unknown  pressure.     As- 
suming isothermal  compression,   the  following  holds: 

PV    =    pv  (5) 

where 
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P  =  unknown  pressure, 

V  =  capture  volume, 

p  =  pressure  in  the  compressed  volume, 

v  =  compressed  volume;    and 

v  =  Ah 


(6) 


wnere 


A    =    capillary  cross -sectional  area,    and 

h     =    difference  in  height  of  mercury  in  the  two  capillaries. 


With 


p    =    P  +  h    (all  pressures  in  mm  Hg), 


Ah 


(7) 


P   =  BX    =    (P-fh)   v 


from  equations  (5)  -  (7). 

In  the  micron  range,   P  «  h  ,    so  that  equation  (8)  becomes 

2 


P    = 


Ah 
V 


=    kh 


(9) 


where    k    =    known  constant. 

In  the  linear  gauge,    a  known  volume  is  captured  and  then 
compressed  to  one  of  three  known  volumes.     The  height  of  a  column 
of  mercury  above  the  compressed  reference  volume  is  then  a  meas- 
ure of  the  original  unknown  pressure.     Assuming  the  isothermal  re- 
lation again,    PV  =  pv  ;  however,   here  v  is  the  known  reference  vol- 
ume; then 

(10) 


v 


P    =    (P+h)^. 


p  =  K  vTv^r) 


(id 
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P    =    Ida  (12) 

where    k   is  a  known  constant. 

CVC  C-TC  100  thermistor  gauges  are  used  to  monitor  the 
stagnation  and  exhaust  chamber  pressures.      These  are  used  to  indi- 
cate pressure  variations  in  the  tank  and  to  check  on  the  autc.inno.ent  of 
steady  conditions,   but  are  not  used  as  a  primary  means  of  pressure 
measurement. 
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HI.     THEORY 
The  flow  of  a  gas  through  orifices,    short  tubes,    and  noz- 
zles at  high  pressure  ratios  in  the  slip  and  transition  regimes  in- 
volves spanning  the  gap  between  gas  dynamics  and  gas  kinetics.     At 
both  ends  of  this  spectrum,    the  governing  equations  and  boundary 
conditions  are  well  defined,    although  solutions  to  the  flow  problem 
may  not  be  readily  obtainable.     For  slip  flow,    the  governing  equations 
are  the  Navier-Stokes  equations,   but  the  boundary  conditions  are  not 
known.     The  conditions  at  the  wall  of  an  object  are  no    longer  the  zero 
velocity  conditions  of  continuum  flow.     In  the  transition  regime,    the 
governing  equations  would  seem  to  be  those  of  gas  kinetics,   but  ac- 
count must  be  taken  of  inter  molecular  as  well  as  molecule -surface 
collisions. 

Continuum  Flow 


Far  into  the  continuum  regime  where  boundary  layers  oc- 
cupy only  a  small  portion  of  the  flow  field,    the  flow  problem  can  be 
analyzed  to  a  good  degree  of  approximation  by  assuming  an'inviscid 
flow  field.     In  this  region,   there  exists  a  definite  critical  pressure 
ratio  (F    /Pa)  beyond  which  the  downstream  conditions  cease  to  in- 
fluence the  upstream  conditions.     The  maximum  mass  flow  is  reached 
at  this  critical  pressure  ratio  for  given  upstream  conditions.      For  a 
smooth  nozzle,   this  condition  is  reached  when  the  Mach  number  is 
unity  at  the  throat,    since  the  maximum  mass  flow  density  occurs  at 
local  sonic  velocity.     This'maxi;-.  pressure  ratio  and  correspond- 

ing mass  flow  are: 
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V 
P  ,    v-T 


m    =    pAV    = 


ly±L 


v     c 


i--%  /x 


2 


V+l 


Although  this  choked  condition  gives  the  maximum  mass  flow,    it  does 
not  determine  the  nature  of  the  flow  downstream  of  the  throat  which 
is  essential  to  optimization  of  a  nozzle.     For  this  analysis,    the 
methods  of  one -dimensional  isentropic  flow  or  characteristics  can  be 
used.     In  these  analyses,   the  boundary  layers  on  the  walls  can  usu- 
ally be  neglected  to  a  good  approximation.     When  they  are  taken  into 
account,    it  is  usually  as  a  small  correction  to  the  nozzle  area  based 
upon  the  boundary-layer  displacement  thickness. 

As  the  slip  regime  is  approached,    the  boundary  layers 
play  a  more  important  role  in  determining  the  over-all  flow.     Even- 
tually they  fill  the  entire  nozzle,    and  the  gas  flow  is  more  like  fully- 
developed  laminar  flow  in  a  long  tube  than  conventional  nozzle  flow. 
Marble  has  done  a  preliminary  analysis  which 

"attempts  to  treat  the  low  Reynolds  number  nozzle  prob- 
lem by  modifying  the  formulation  to  include  the  viscous 
stresses  associated  with  wall  shear,   but  neglecting  the 
normal  viscous  stress  associated  with  acceleration 
along  the  flow  direction.     It  is  assumed  that  the  flow  is 
dissipative,   but  locally  adiabatic.     Under  these  restric- 
tions the  results  should  be  adequate  to  indicate  the  order 
of  magnitude  of  the  losses  and  of  their  effects  on  nozzle 
geometry.  "  2 
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In  this  analysis,    the  velocity  dependence  on  axial  distance 
along  the  nozzle  is  chosen.      This,    together  with  the  gas  thermody- 
namic properties  and  an  estimate  of  the  cross-sectional  velocity 
profile, yields  the  information  required  to  estimate  the  nozzle  ge- 
ometry to  achieve  a  given  nozzle  performance.     Marble  concludes 
that  viscous  dissipation  plays  a  strong  role  for  small  Reynolds  num- 
bers,   and  that  for  a  given  mass  flow,    as  the  length  of  the  nozzle  is 
increased,    the  area  ratio  must  increase  to  achieve  the  same  dis- 
charge velocity.     For  very  low  Reynolds  numbers,    the  nozzle  di- 
vergence angle  must  be  increased  as  the  length  is  increased. 

Although  the  choked  condition  for  a  nozzle  in  continuum 

flow  occurs  at  modest  pressure  ratios  (P    /P     =  1.  89  for  air  and 

r  c      e 

2.  06  for  He),    much  higher  ratios  are  required  for  the  flow  through 
an  orifice.     This  is  because  the  sonic  line  for  an  orifice  is  S-shaped, 
while  that  of  the  nozzle  is  planar.     Since  the  flow  in  the  plane  of  the 
orifice  is  not  sonic,    it  follows  that  the  mass  flow  is  smaller  than 
for  a  smooth  nozzle  of  the  same  throat  area. 

The  two-dimensional  analogy  to  orifice  flow,    flow  through 

3 

a  slit,    has  been  analyzed  by  the  method  of  characteristics    .     The 

flow  becomes  sonic  at  the  lip  of  the  slit  and  flows  locally  around  the 
lip  like  a  Prandtl- Meyer  expansion.     The  streamline  maps  on  a 
characteristic  in  the  hodographic  plane.     When  the  characteristic 
connects  the  sonic  point  on  the  axis  with  the  corner,    the  flow  down- 
stream has  no  effect  on  the  upstream  conditions  and  the  flow  is 
choked.     This  characteristic  is  an  epicycloid  from  the  sonic  circle  on 
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the  u-axis  and  meets  the  corner  characteristic  at   8  =  45   .      The 
critical  pressure  ratio  is  the  one  corresponding  to  this  45     turn  from 
sonic  in  Prandtl- Meyer  flow.     Liepmann  gives  the  critical  pressure 
ratios  as  58.  3,    25.  6,    and  21.  5  for  y's  of  5/3,    7/5,    and  4/3,    re- 
spectively . 

Obtaining  the  mass  flow  is  more  difficult,    for  the  full  hodo- 
graph  equation  must  be  solved.     Frankl  has    computed  for  a  two- 
dimensional  slit  that  the  ratio  of  the  mass  flow  to  the  mass  flow  of  a 

5 
smooth  nozzle  of  equal  area  is  0.  85  for  a  gas  with  y  =  1.  4. 

The  axisymmetric  problem  of  the  orifice  is  even  more  dif- 
ficult since  the  hodograph  equations  are  not  linear.     Liepmann  gives 
an  argument  that  shows  that  the  required  Prandtl-Meyer  turning 
angle  and  pressure  ratio  are  larger  for  an  orifice  while  the  mass 
flow  is  smaller  than  for  a  slit.     The  differences  between  the  two 
should  be  small.     Thus,    it  is  expected  that  the  mass  flow  through  an 
orifice  at  high  Reynolds  number  will  be: 

1    y+l 


m  =  a-4g--\/£ 


*     c 

4 


12  Y-l 


Y+l 


A    ,  (15) 


where    a  is  on  the  order  of  0.  85. 

As  the  Reynolds  number  is  decreased,   boundary  layers 
build  up  near  the  lip  and  into  the  orifice.     These  can  have  two  effects. 
First,    they  reduce  the  effective  area  of  the  orifice,    and  second,    they 

can  round  out  the  orifice  entrance  and  tend  to  increase  the  mass  flow. 

4 
Liepmann     states  that  a  limiting  formula  at  high  Reynolds  number 

should  have  the  form 
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i 


m    =    mD    ^      (1  -  CRe    2)  (16) 

Re  >  oo  x      ' 

where   C   is  a  constant.     He  states  that  experimentally   C   is  of  the 
order  of  \  . 

The  third  profile  of  interest,    the  short  tube,   would  be  ex- 
pected to  act  similarly  to  an  orifice  in  the  continuum  regime  with  the 
same  entrance  and  boundary  layer  effects. 

Free-Molecule  Flow 

In  the  free -molecule  regime,    where  the  mean  free  path  is 
much  larger  than  the  characteristic  dimension  of  the  profile  in 
question,   the  flow  through  the  profile  can  be  determined  by  the  meth- 
ods of  kinetic  theory.     It  is  a  matter  of  statistically  keeping  track  of 
molecules  and  their  surface  interactions. 

Following  Present   ,    one  can  compute  the  number  of  mole- 
cules passing  through  an  orifice.     Consider  an  element  of  area  on  the 
wall  dA  .     Choose  a  system  of  spherical  coordinates  with  origin  at 
the  center  of  dA  and  axis  normal  to  dA  .     The  mean  free  path  of  the 
gas  molecules  is 

KT 

z r 

,2      =     -,/T 

rra 


x   =      7=    "     2    =      r     z  '  (17) 

"V2P  Trd^        V2 


The  mean  molecular  velocity  is 

c    =    V(8/tt)  RT  .  (18) 


c 

The  collision  frequency  per  molecule  is 

2 


F    h.  V2  Tracn  (19) 
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where   n  =  molecule  number  density  and  a  =  collision  diameter.     The 
number  of  molecule  collisions  in  a  volume  element  dv  located  at   r  , 
0  ,    and   <{>  is  equal  to  the  average  number  of  molecules  in  dv  times 
the  probability  per  unit  time  that  a  molecule  collides, 

ndvF    =    (c/\)ndv    .  (20) 

Molecules  leave   dv  in  random  directions.     The  fraction  that  leaves 
dv  in  the  direction  of  dA  is    dA|  cos  9  |  /4r     .     The  fraction  of  mole- 
cules leaving  dv  toward  dA  which  arrive  at  dA  without  undergoing 
further  collisions  is 

cndy     dA|  cos  8  1      -r/\ 
*  4„2  C 

-r/X 
where   e  is  the  probability  that  a  molecule  travels  a  distance    r  , 

from  its  last  collision,   without  a  further  collision.     The  total  number 

of  molecules  striking  a  unit  area  of  wall  per  unit  time  is 


N    =     cn 


4tT\         y 


I 


-^ele"r/Xdv  (21) 


cn 
4trX    , 


r 

2ir        tt/2  co 

d<j>    f  d9  sin9  cosG     |dre"r'X 


Now,    consider  a  small  hole  of  diameter   d   much  smaller 
than  the  mean  free  path.     A  molecule  in  the  vicinity  of  the  hole  and 
traveling  toward  it  will  pass  through  the  hole  without  suffering  any 
collision  in  the  neighborhood.     The  molecules  pass  through  individu- 
ally  and  independent  of  the  motion  of  other  molecules   • 
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Equation  (22)  gives  the  number  flow  of  molecules  through 
a  sharp -edged  orifice  when  the  Knudsen  number  is  sufficiently  high. 
It  can  be  applied  to  the  gas  on  both  sides  of  the  orifice  to  give  the  net 
flow,    since  the  two  flows  are  independent  of  each  other  if  the  mean 
free  paths  on  each  side  of  the  orifice  are  much  larger  than  the  diam-- 
eter.     This  type  of  flow  is  called  effusive  flow  or  effusion.     If  the  gas 
on  each  side  of  the  orifice  is  at  the  same  temperature,   the  net  mass 
flow  rate  is  therefore 

1 


m    = 


V2ttRT 


(Pc-Pe)  <23) 


where 


A    =    area  of  the  orifice, 
P      =    upstream  pressure, 

P       =    downstream  pressure. 

e  r 

When  the  orifice  is  not  sharp -edged  but  has  a  finite  length, 
some  of  the  molecules  collide  with  the  side  walls  of  the  orifice  and 
then  are  re -radiated  in  a  random  fashion.     This  causes  some  of  the 
molecules  that  enter  the  orifice,    or  short  tube,   to  be  re-emitted  to 
their  original  side  of  the  orifice.     Clausing  showed  that  the  correct 
expression  for  circular  cross-sections  is 

rh    =  * A(Pc-Pe)  (24) 

"V2ttRT  c      e 

where    x  is  a  function  of  the  length-to-diameter  ratio  of  the  orifice 

or  tube  and  is  called  the  Clausing  factor.     As  the  tube  becomes  long, 

7 
H  approaches    4d/3«t  as  an  upper  limit.     Present  gives  a  derivation 

for  the  mass  flow  for  long  tubes  and  shows  that   h=  4d/3£'from  geo- 
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metric  and  gas  kinetic  considerations    . 

For  long  tubes  of  varying  cross  section  A  and  perimeter   H  , 

the  fundamental  relation  deduced  by  Knudsen  is  (cs.    refs.    7,  8) 

•     _    4           c                   c       e  .-,£., 

m    "    3    T — RT *  {Z5) 


B 


dx 

.        oA/ 

For  a  conical  shape  of  minimum  diameter   d  ,    half  angle    8  ,    and  axi- 
al length  t  $    this  relation  gives 

I  2 

4  cd      v>    w  8*  w 1      y  d    (d/2  +  £tan9)  /9,x 

m    =    3  (Pc"Pe)(RT)(-2T7d)2        (d  +  <,taa6).        •  {26) 

The  flow  for  other  shapes  in  the  free-molecule  regime  can 

likewise  be  calculated  from  geometric  and  gas  kinetic  considerations 

(cs.    refs.    9»  10). 

Slip  and  Transition  Flow 

For  flow  in  the  slip  and  transition  regimes,    there  exist  a  few 
formulae  to  calculate  the  flow  through  orifices  and  tubes.     These  are 
in  the  form  of  empirical  or  semi-empirical  corrections  to  the  con- 
tinuum equations.     As  an  example,    consider  Poiseuille  flow  of  a  gas 
through  a  long,    straight  pipe  under  the  restrictions  that  (a)  the  gas  is 
incompressible;  (b)  the  flow  is  fully  developed,    i.  e.  ,    the  cross- 
sectional  velocity  profile  is  constant  throughout  the  length;  (c)  there 
is  no  turbulent  motion;  and  (d)  the  flow  velocity  at  the  tube  wall  is 
zero.     The  gas  effectively  flows  in  concentric  cylinders  which  slide 
over  each  other.     Since  the  flow  is  steady,    the  viscous  shearing  drag 
on  the  curved  surface  is  balanced  by  the  pressure.     Therefore, 
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2irrdxu  ~    =    -rrr2dP  .  (27) 

dr 


With  the  velocity  zero  at  the  walls,  integration  of  equation  (27)  gives 

(28) 


u 


d2-r2     dP 


16  |j,         dx 

The  total  molecular  flow  (N)  through  the  tube  is 

d/2 

N    =    n 


*  ^4 


o 


u2irrdr    "    T2^lKT  P  "dx"    • 


and 


N    = 


,4  ,P  2-P  2 

rrd  /     c        e 


128^iKT<C 


(-Vs-)  <29> 


since   N  is  independent  ofx  by  conservation  of  molecules.    The  mass 

flow  rate  is  therefore 

4   -  [ 

TTd^   P     AP  ,_n. 

m    "    728  U  RTT  (30)  • 

where   P  =  (P  +P    )/2   and  AP  =  P    -P     . 
c      e  c      e 

Although  this  Poiseuille  equation  is  strictly  limited  by  re- 
strictions given  above,    it  has  been  modified  for  use  in  slip  flow.     The 
argument  for  this  is  given  below.     However,    since  the  correction  is 
based  on  experimental  data,    it  should  be  pointed  out  that  the  modified 
equation  must  include  a  correction  for  compressibility  effects,    which 
are  present,    even  though  the  argument  is  based  solely  on  slip. 

When  the  mean  free  path  is  smaller  than  the  characteristic  di- 
mension,  but  not  negligibly  small,    an  additive  term  to  represent  the 
slip  at  the  surface  is  used.     It  is  assumed  that  the  gas  at  the  wall  has 

a  finite  velocity  u     .     The  number  of  molecules  striking  a  unit  area  of 
7      o  ° 

wall  surface  per  unit  time  is 


N    =    n 
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+  /     d      \  dP 
uo      \  128|i  /  dx 


(3D 


An  estimate  is  made  for  u     from  mean  free  path  considerations.     It 
is  found  to  be  . 

2\    /du\  Z\     d    dp  ,„-,» 

u      =    — r-  It-J  =    -5-  T~  "^n*  (32; 

o  3     Ndr'  3    4u   dx 

r=a  K 

where    r  =  cylindrical  radial  coordinate  from  the  axis  of  the  tube,    and 

x  =  cylindrical  coordinate  along  the  tube.     The  modified  Poiseuille 

flow  equation  is  (cf.    ref.    6) 

4  3 

•     -         ffd  ^  AP        ud       AP  ,,,. 

m    "    T2^lRT    ^~T  +  ~8c~  T"    '    .  l      ' 

The  additive  term  has  the  same  form  as  the  free-molecule  mass  flow 
through  a  tube,   but  is  smaller  by  a  factor   3tt/16.     Equation  (33)  is 
valid  over  a  limited  range  and  fails  when  the  mean  free  path  and 
characteristic  dimension  are  approximately  equal. 

■ ,  Dushman  gives  slip  flow  formulae  for  tubes  attributed  to  ' 
three  sources.     From  Kennard  he  gives 


rrd4      «■   AP 


m   =      "r   .    P 


128  (j£  RT 


(1+i£)  <34> 


where   Q    is  the  coefficient  of  slip  given  by 

£=¥^(-^)  '       <«> 

and  f  is  the  transfer  ratio  of  momentum  and  represents  the  fraction 
of  molecular  collisions  with  the  walls  which  result  in  diffuse  scat- 
tering of  the  molecules.     The  value  of  f  varies  between  about  0.  8 
and  1.  0  . 
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Knud8en  found  that  the  flow  through  a  long  tube  can  be 

written 

rh    =    rh     +  Zrh.  (36) 

c  fm 

where   rh      is  the  continuum  flow  rate  and   mr         is  the  free  molecule 
c  fm 

value.     Z   is  a  function  of  the  mean  pressure  in  the  tube,   the  tem- 
perature,   the  tube  radius,    and  the  viscosity  of  the  gas. 

i  +  i.z35iyo^P 


Sreekanth  has  recently  developed  a  modified  Poiseuille 
equation  from  experimental  data  on  the  flow  through  short  tubes  (-6/d 
from  0.  005  to  0.  995)  at  pressure  ratios  from  approximately  1:1  to 

30:1  and  at  Knudsen  numbers  from  0.  1  to  2  (the  free -molecule  end  of 

12 
the  transition  regime)      .     This  modified  equation  is 


m    = 


^4       gAP  +  0-519d3  ^ 


128(jRT 


1 


which  is  almost  identical  to  the  Poiseuille  slip  flow  relation,    equation 
(33),    divided  by  the  factor  (1  +  d/l)  . 

The  above  equations  illustrate  some  of  the  empirical  rela- 
tions that  cover  flows  in  parts  of  the  slip  and  transition  regimes.     No 
known  empirical  relation  or  theory  covers  the  entire  range  from 
free -molecule  to  continuum  flow.- 
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IV.     RESULTS 

Experiments  were  conducted  to  determine  the  mass  flow 
through  an  orifice  with  an  £/d   of  0.  0895  (Figure  9)  at  Knudsen 
numbers  of  from  8.  2  to  0.  00068.      The  pressure  ratios  across  the 
orifice  were  at  least  100:1  .     The  experiments  were  conducted  to 
verify  that  the  low-density  gas  dynamic  facility  was  functioning  prop- 
erly and  that  the  range  from  continuum  to  free -molecule  flow  could 
be  covered. 

The  initial  set  of  data,   which  took  one  month  to  collect, 
revealed  that  there  was  a  slight  leak  into  the  stagnation  chamber. 
The  leak  rate  was  approximately  one  micron-liter  per  second.     This 
caused  the  measured  flow  rates  to  be  too  low  near  and  in  the  free- 
molecule  regime.     The  measured  flow  rates  dropped  to  about  85  to 
90  per  cent  of  the  expected  value.     This  small  leak  had  virtually  no 
effect  on  the  data  collected  in  the  continuum  regime. 

The  initial  test  also  revealed  some  difficulties  with  the 
flow  measurement  system  and  over-all  system  stabilization.     When 
the  Vol- U- Meter  was  used  for  flow  measurement,    the  1  per  cent 
pressure  change  in  the  input  to  the  variable  leak  when  the  flow  was 
passed  through  the  Vol-U-meter  caused  an  upset  in  the  steady-state 
previously  attained.     This  was  remedied  by  throttling  the  gas  in  the 
feed  line  in  parallel  with  the  Vol-U-Meter  to  match  the  expected 
output  pressure  of  the  Vol-U-Meter.     This  meant  that  no  change  in 
feed  pressure  occurred  when  the  flow  was  directed  through  the  Vol- 
U-Meter. 
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At  high  Knudsen  numbers  it  was  found  that  it  takes  as 
much  as  five  hours  for  the  system  to  reach  steady-state.      This  is 
due  to  the  relatively  large  volumes  and  small  flow  rates.     After  the 
system  stabilizes,    it  takes  about  one  hour  to  collect  one  data  point. 
The  system  is  sensitive  to  temperature  changes,    the  maximum  al- 
lowable rate  being  about  2   F/hr. 

As  the  flow  rate  increases  with  decreasing  Knudsen  num- 
ber,   the  time  constant  of  the  system  decreases  and  steady-state  is 
attained  fairly  rapidly  --  in  as  little  as  fifteen  minutes. 

After  the  initial  difficulties  were  remedied,    another  series 
of  data  collecting  runs  were  made.     A  typical  test  run  is  outlined  in 
Appendix  I.     It  is  these  data  that  will  be  discussed  in  the  following 
sections. 

Data  Reduction 

Two  methods  were  used  to  reduce  the  data  so  that  the 

results  could  be  compared  with  the  results  of  Liepmann     and 

12 
Sreekanth 

Liepmann  defines  a  characteristic  velocity,    W  ;  Mach 
number,    Ma  ;  and  Reynolds  number,    Re  ;  by 

P   -P     * 


W    "    (-T-*)  <39> 


Pc 

P    -P     - 

**  -  HS  =  (-v*)*  wo) 

(Pc/P>  *c 

Re   =  .24   =  -4.  (lE-ll)1  (41) 

c  c  rc 
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where   v      is  the  kinematic  viscosity  of  the  gas  evaluated  at  the  up- 
stream conditions.      He  writes  the  mass  flow  rate  in  terms  of  the 
orifice  area,    A  ,    and  characteristic  velocity,    W  ,    in  the  form 

m    =    Tp   WA  (42) 

where  T   is  a  dimensionless  function  which  depends  on: 

"Ma  ,    Re  ,    the  ratio  of  specific  heats    y  ,    and  in 
the  transition  region  possibly  on  slip,    accomoda- 
tion coefficients,    and  possibly  a  Reynolds  number 
based  on  the  bulk  viscosity  coefficient.  "4 

The  parameter  T   is  well  defined  in  the  two  limits  of  con- 
tinuum and  free-molecule  flow.     The  free-molecule  limit  from 
Knudsen's  work  is: 

rK  =  -^-  . '  (43) 

K      VI7 


In  the  continuum  limit  with   Ma  =1.0, 


r    =    a(y) 


2+1' 

y-1 


LVW+: 


(44) 


where   a(y)   is  theoretically  well  defined, varying  from  a  =  1   for  a 
smooth  nozzle  to  approximately  a  =  0.  85  for  an  orifice. 

When  the  Clausing  factor  (k)  is  included  to  modify  the 
free-molecule  limit,    the  relation  for  T  becomes 

T*-it     ■ 

JLiepmann  presents  his  results  in  the  form  T/TK  versus    I /Re  *     The 
continuum  limit  for  the  former  function  is 
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Sreekanth  presents  his  results  in  the  form  rh/rh, 

fm 

versus  Knudsen  number  defined  as  the  ratio  of  mean  free  path  evalu- 
ated at  stagnation  chamber  conditions  to  diameter  of  the  orifice  or 

12 

short  tube.      The  expression  used  to  calculate  mean  free  path  is 

c 
As  given  by  equation  (24),    the  free -molecule  mass  flow  limit  is 


i— XA(P   -pJh    . 


m,. 

lm  (2TrRTc)t 

This  means  that  the  continuum  limit  for  the  dependent  variable  is 

1    V+l 
2  " 


m  i 

— s.  ■  =  a  <2»V)* 

mfm 


V+l 


2  v-1  x 

(1  -  P    /P    )      '  *  (48) 


In  the  limit  with  P     >»  P     ,    this  reduces  to  Liepmann's  T jTy,  . 

c  e  is. 

It  should  be  noted  that  the  independent  variables  of  Liep- 

mann  and  Sreekanth  are  closely  related,    and  in  fact  differ  by  a  small 

constant  factor  in  the  limit  of  P     >»  P     .     The  ratio  Kn  to    1/Re   is: 

c  e 

v  16u  RT      k  .  P   -P     %  ',  P 

This  reduces  to  1.  275  when  P     >»  P     . 

c  e 

Thus,    the  two  systems  of  data  reduction  are  similar,   but 
not  exactly  the  same.     Both  approach  the  same  theoretical  limits  for 
continuum  and  free-molecule  flow.     The  present  data  were  reduced 
by  both  of  these  schemes.. 

Data  Presentation  and  Analysis 

Figure  11  shows  the  present  data  superimposed  on  Liiep- 
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mann's.     The  two  sets  conform  above    1/Re  =10     '  and  both  seem  to 
be  slowly  approaching  the  same  free-molecule  limit  of  1.  0  .     At  the 
lower    1/Re   end,   the  two  sets  seem  to  be  asymptotically  approaching 
different  limits.      This  is  as  would  be  expected.      The  present  data 
were  for  an  orifice  with  <t/d  =  0.  0895,   while  I_.iepmann's  orifice  had 

<t/d  =  0.  0254.     The  Clausing  factors  for  the  two  are,    respectively, 

7 
0.  918   and   0.  975.        Assuming  that  a  =  0.  85  ,    the  two  theoretical 

values  ibr    T/TK  at  the  continuum  limit  are  1.  68  and  1.  56,    resocc-    • 

tively.     The  last  few  values  of  T/Pt,   obtained  from  the  present  data 

were  equal  to  the  theoretical  limit  of  1.  68  .     Thus,   both  sets  of  data 

appear  to  be  compatible  at  both  limits  and  throughout  the  slip  and 

transition  regimes.     They  both  show  the  same  smooth  transition 

from  continuum  to  free -molecule  flow. 

-2  -1 

Between    1/Re  =  10    "   and    10       ,    lAepmann's  data  exhibit 

a  maximum  which  he  attributes  to  the  boundary  layer  smoothing  the 
entrance  condition  to  the  orifice  and  thereby  increasing  the  mass 
flow.     The  present  data  do  not  exhibit  this  maximum,   but  do  exhibit 
a  local  maximum  or  a  change  in  curvature  which  can  be  explained 
qualitatively  by  the  same  boundary  layer  argument. 

Figure  12  shows  the  present  data  together  with  Sreekanth's 
data  for  an  orifice.     The  similarity  between  this  presentation  and 
Figure  11  is  evident.     Again,   both  sets  of  data  are  compatible. 
Sreekanth's  data  do  not  extend  far  enough  into  the  low  Knudsen  num- 
ber range  to  get  a  good  comparison  throughout  the  slip  and  transition 
ranges,   but  a  comparison  at  the  free-molecule  end  shows  good  agree- 
ment. 
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Sreekanth's  empirical  formula,    equation  (38),    is  also 
shown  in  Figure  12.     It  is  satisfactory  at  the  high  Knudsen  number 
end  of  the  transition  regime.     It  overestimates  the  flow  rate  below  a 
Knudsen  number  of  approximately  0.  2  . 

Over  all,    the  present  data  show  the  transition  from  continu- 
um to  free-molecule  flow.     The  majority  of  the  transition  takes  place 
within  two  orders  of  magnitude  (0.  01  <  Kn  <  1  ,    or    100  <  Re  <  1   ). 
Boundary  layer  and  slip  effects  seem  to  be  evident  in  the  range 
0.  001  <  Kn  <  0.  01  ,    or    1000  <  Re  <  100  .     The  slow  approach  to  the 
free-molecule  limit  extends  beyond  Kn  =  10   or  Re  =  0.  1  . 
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V.     CONCLUSIONS 

The  continuous  flow,    low-density  gas  dynamic  facility  de- 
signed to  study  gas  flow  through  various  profile  apertures  in  the 
slip  and  transition  regimes  was  successfully  designed,    constructed, 
and  tested.      The  test  was  performed  with  an  orifice  as  the  test  pro- 
file.     The  data,    which  exhibit  surprisingly  little  scatter,    agree  with 
previous  data  by  Liepmann  and  Sreekanth.      The  facility  can  be  used 
to  gather  data  over  a  larger  range  than  any  previously  reported. 
The  range  extends  from  the  threshold  of  free-molecule  flow  to  con- 
tinuum flow. 

The  collected  data  show  that  the  transition  to  continuum 
from  free-molecule  flow  takes  place  in  essentially  two  orders  of 
magnitude  in  Knudsen  number.     The  transition  is  smooth,   with  no 
discontinuities  as  would  be  associated  with  an  abrupt  change  in  the 
flow  field.     A  comparison  of  the  data  with  the  most  recent  empirical 
formulae  shows  that  the  re  is  no  universal  equation  to  cover  the  entire 
spectrum.     The  measurements  show  that  the  flow  is  approaching  the 
theoretical  limits  at  both  high  and  low  Knudsen  numbers. 

Some  interesting  questions  were  raised  by  this  study.     The 
region  around   Kn  =  0.  01    or   Re  =  100   exhibits  a  local  maximum  or 
change  in  curvature  that  Liepmann  reasoned  was  caused  by  boundary 
layers  smoothing  the  orifice  entrance  conditions.     This  area  is 
worthy  of  a  more  detailed  investigation.     The  transition  to  free- 
molecule  flow  does  not  appear  to  be  complete  at  a  Knudsen  number  of 
10.     The  final  approach  to  the  theoretical  free-moiecule  limit  appears 
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i 

to  occur  very  slowly.     The  final  attainment  of  this  theoretical  limit 
also  warrants  further  study. 

In  conclusion,    it  can  be  said  that  the  facility  meets  its 
design  criteria.     It  is  hoped  that  it  will  be  used  to  investigate  the 
areas  of  orifice  flow  mentioned  above  and  flows  through  other  pro- 
files such  as  short  tubes  and  nozzles. 
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APPENDIX  I.     Typical  Test  Procedure 
I.     Initial  Pump  down 

A.  Commencing  with  all  valves  closed: 

1.      Open  shutoff  valves  to  vacuum  pressure  measurement 

lines  (3  valves) 
2..    Open  tank  bypass  valve 
3.     Open  4"  Temescal  gate  valve 

B.  Start  cooling  water  flow  to: 

1.  Stokes  pump  --  approximately  2  gal/min 

2.  CVC  diffusion  pump  --  approximately  1/15  gal/min 

C.  -  Start  all  pumps: 

1.  Stokes  pump 

a.  Check  solenoid  oil  valve  with  ferromagnetic  ma- 

terial 

b.  If  not  magnetized,    shut  down  immediately 

2.  Cenco  Hypervac  25  backing  pump 

3.  CVC  diffusion  pump 

D.  Fill  McLeod  gauge  cold- trap  reservoir  with  liquid  nitrogen 

E.  Open  McLeod  gauge  stopcocks 

F.'     Close  variable  leak  to  a  counter  setting  of  010. 
G.     Open  shutoff  valve  downstream  of  variable  leak 
H.     When  exhaust  chamber  pressure  as  read  by  thermistor  no.    1 
is  10  microns: 

1.  Close  4"  Temescal  gate  valve 

2.  Open  2"  Temescal  gate  valve 
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3.     Shut  down  Stokes  pump 
I.      Pump  down  system  to  desired  ultimate  pressure.       Note:    if 
this  is  the  first  pumpdown  after  the  system  has  been  ex- 
posed to  atmospheric  conditions: 

1.  Open  Stokes  pump  gas  ballast  valves  (2) 

2.  Delete  steps  D  and  E 

3.  Operate  on  gas  ballast  for  15  -  20  minutes 

4.  Close  gas  ballast  before  going  to  Step  H 

5.  Pump  down  system  for  at  least  24  hours  before  pro- 

ceeding (it  may  be  desirable  to  pump  down  for  a  few 
days  to  minimize  outgassing) 

6.  Execute  steps  D  and  E  before  proceeding 
II.     Test  Gas  Flow  Setup 

.A.     Close  tank  bypass  valve 

B.  Close  one  of  the  shutoff  valves  to  the  McJLeod  gauge  vacuum 
-line  (which  one  depends  on  whether  exhaust  or  stagnation 

chamber  pressure  is  to  be  measured) 

C.  Fill  feed  line  cold-trap  reservoir  with  liquid  nitrogen 

D.  Start  constant-temperature  water  bath 

1.  Start  motor  and  heater 

2.  Set  temperature  regulator  after  flow  rate  is  established 

E.  Open  all  valves  in  the  gas  feed  line  except: 

1.  Close  individual  rotameter  shutoff  valv 

2.  Close  rotameter  cutlet  valve 

3.  Close  feed  line  shutoff  valve  immediately  downstream 

of  the  £eed  line  vent  valve 
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F.  Open  helium  (gas  supply)  bottle 

1.  Set  high-pressure  regulator  to  about  100  pci  (must  be 

below  200  psi) 

2.  Set  low-pressure  regulator  to  about  5"  K_0 

G.  After  a  few  seconds  of  venting  feed  gas  to  the  atmosphere: 

1.  Close  feed-line  vent  valve 

2.  Open  feed-line  shutoff  valve  dov/nstream  of  the  vent 

H.     Open  variable  leak  to  the  desired  counter  reading  (graph  of 
approximate  counter  reading  versus  stagnation  pressure  is 
.  helpful  for  this ) 

1.  Monitor  exhaust  chamber  pressure  on  thermistor  no.    1 

2.  If  exhaust  chamber  pressure  exceeds  10  microns, 

switch  pumps 

a.  Close  2"  Temescal  gate  valve 

b.  Start  Stokes  pump  and  wait  30  seconds 

c.  Open  4"  Temescal  gate  valve 

I.       When  satisfied  with  pump  setup,    monitor  stagnation  cham- 
ber pressure  on  thermistor  no.    2,    0-50  mm  dial  gauge, 
or  manometer.     Note:    Keep  manometer  shutoff  valve 
closed  below  50  mm  Hg  to  minimize  mercury  vapor  con- 
tamination of  the  variable  leak. 
III.     Flow  Rate  Measurement 

A.     Low  Flow  Rate  --  Vol -U- Meter 

1.  Open  input  valve  to  rotameter  no,    1 

2.  Open  rotameter  output-line  shutoff  valve 
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3.  Close  rotameter  bypass  valve 

4.  Partially  open  rotameter  no.    1  needle  valve  to  give 

about  .  1  psi  pressure  drop  between  flow  measure- 
ment system  input  pressure  (0-15  psia  dial  gauge) 
and  Vol-U- Meter  output  pressure  (0  -  25,    25  -  50  psia 
dial  gauge) 

5.  Stabilize  system  for  at  least  2  hours 

6.  Continually  refill  cold  traps 

7.  Close  Vol-U- Meter  bypass  valve  and  simultaneously 

open  rotameter  bypass  valve.     Note:    this  procedure 
minimizes  changes  in  the  gas  input  pressure  to  the 
.  •   variable  leak  when  the  Vol-U- Meter  is  thrown  into  the 
system. 

8.  Time  the  piston  between  any  two  volume  marks  (pref- 

erably 0  to  25  cc  at  the  slower  rates,    and  0  and  22  cc 
at  the  faster  rates) 

9.  Before  the  piston  reaches  the  top  of  the  tube: 

a.  Stop  the  piston  by  cracking  open  the  Vol-U- Meter 

bypass  valve 

b.  Lower  the  piston  gently  to  the  bottom  of  the  tube  by 

carefully  opening  the  Vol -U- Meter  bypass  valve  a 
bit  more 
Bt     Moderate  to  High  Flow  Rates  --  Rotameters 

1.  Close  input  needle  valve  of  desired  rotameter 

2.  Open  input  valve  to  rotameter 
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3.  Open  rotameter  output-line  shutoff  valve 

4.  Close  rotameter  bypass  valve 

5.  Open  rotameter  input  needle  valve  carefully  to  full 

open  position.     Note :     this  should  be  done  carefully  to 
prevent  damage  to  the  rotameter  float. 

6.  The  rotameter  can  be  read  continuously  and  will  help 

indicate  when  the  flow  is  stabilized  (the  primary  means 
of  assuring  stabilization  is  stagnation  chamber  pres- 
sure measurements). 
IV.     Other  Data  Measurements 

A.  Atmospheric  conditions 

1.  Pressure  -  barometer  (apply  temperature /gravity  cor- 

rection to  indicated  Hg  column  height) 

2.  Temperature  -  thermometer  on  barometer  case 

B.  Gas  inlet  conditions 

li     Pressure  -  before  flow  measurement  system 

a.  Low-pressure  regulator  dial  gauge  (0  -  30  in.  H-O) 

-  this  pressure  plus  atmospheric  should  equal  flow 
measurement  system  input  pressure 

b.  Flow  measurement  system  input  dial  gauge  (0-15 

psia) 

c.  Flow  measurement  system  input  manometer  (0  - 

800  mm  Hg  gauge) 

3.  Temperature  -  before  flow  measurement  system 
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C.  Stagnation  Chamber  Conditions 

1.  Pressure 

a.  Thermistor  no.    2 

b.  McL.eod  gauges 

c.  Dial  gauge  (0  -  50  mm  Hg) 

d.  Manometer  (0  -  800  mm  Hg)  -  the  scale  is  already 

corrected  for  temperature  and  gravity  -  subtract 
scale  reading  from  corrected  barometer  reading 

2.  Temperature  -  thermometer  is  in  protective  shroud 

protruding  from  upstream  end  of  tank 

D.  Exhaust  Chamber  Conditions  -  Pressure 

1.  Thermistor  no.    1 

2.  McJLeod  gauges 
V.     McL»eod  Gauge  Operation 

A.  Activate  pressurization  system 

1.  Close  individual  gauge  shutoff  valves 

2.  Open  vent  line  valve 

3.  Open  needle  control  valve 

4.  Open  nitrogen  bottle  valve 

5.  Open  regulator  shutoff  valve 

6.  Set  regulator  pressure  at  approximately  25  psig 

7.  Close  needle  control  valve 

B.  Sample  measurement 

1.     Capture  of  gas  sample  (either  gau  ;e) 

-.  a.     Close  vent  line  valve 
b.     Open  desired,  gauge  shuccff  valve  1  3/4  turns 
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c.  .  Open  needle  valve  slor/ly  to  obtain  rate  of  mercury 

rise  such  that  sample  capture  takes  15  -  30  sec- 
onds 

2.  Non-linear  gauge 

a.  After  sample  capture,    increase  rate  of  mercury 

rise  but  slow  it  drastically  near  the  top  of  the 
large  volume  to  avoid  shocking  the  closed  capil- 
lary tube 

b.  When  mercury  is  in  the  closed  capillary,    increase 

the  rate  of  rise  again 
C.     Slow  the  mercury  rise  as  the  level  in  the  open  capil- 
lary nears  the  top 

d.  Stop  the  top  of  the  mercury  meniscus  at  the  bottom 

of  the  black  tape  line  by  either  closing  the  needle 
valve  or  closing  the  gauge  shutoff  valve 

e.  Tap  both  capillaries  a  few  times  to  overcome  sur- 

face tension  effects 

f.  Read  the  difference  in  mercury  column  heights 

3.  Linear  gauge 

a.  •  After  sample  capture,    increase  rate  of  mercury 

rise  but  slow  it  near  the  top  of  each  large  volume 
to  avoid  overshooting  the  compression  volume  cal- 
ibration mark 

b.  Stop  the  mercury  meniscus  a-  the  desired  scribed 

volume  mark  by  closing  ci:/.    r  :he  ^auge  shutoff 
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vaAvt*  or  the  needle  valve 

(1)  Use  1st  scribe  mark  and  left  tube  for  pres- 
sures between  125  and  10  mm 

(2)  Use  2nd  scribe  mark  and  middle  tube  for 
pressures  between  10  and  1  mm 

(3)  Use  3rd  scribe  mark  and  right  tube  for  pres- 
sures .between  1  and  .  1  mm  (pressures  below 
.  15  mm  can  also  be  read  on  the  non-linear 
McLeod) 

c.     Read  the  difference  in  mercury  column  ho  i^  .vis 

4.  To  lower  the  mercury  level  in  either  gauge: 

a.  Close  needle  valve 

b.  Open  vent  line  valve 

c.  Open  gauge  shutoff  valve 

Mote:    if  a  slug  of  mercury  hangs  up  in  the  closed  capil- 
lary of  the  non-linear  gauge,    carefully  heat  the  tube 
with  a  fuel-rich  (yellow)  natural  gas /oxygen  flame  until 
the  slug  vaporizes. 

5.  Conversion  from  scale  reading  to  pressure 

a.  Apply  temperature  and  gravity  correction  to  all 

mercury  column  heights  before  using  scale  factors 

b.  .   Scale  factors 

(1)  0  -  100  micron  gauge:    P.    x  =  .  2375XlO_Zh       2 
•    '    y  B     *  (|i)  mm 

-3 

(2)  0-1  mm  gauge:    P.         •  =-l.LL,'>~iZ      h 

°      °  kmm)  mm 

(3)  0-10  mm  gauge:    P.         .  =  1.894X10     h 

■  °      °  .nv.-n,/  mm 
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(4)    0-  100  mm  gauge:   P,         ,  =  2.  494X10"  1h 

fa     °  (mm)  mm 

VI.     System  Shutdown 

A.  Test  gas  feed 

1.  Close  valve  downstream  of  variable  leak 

2.  Close  feed-line  shutoff  valve  downstream  of  feed  line 

vent 

3.  Check  Vol- U- Meter  bypass  and  rotameter  bypass 

valves  open 

4.  Turn  off  helium  (test  gas)  bottle 

5.  Open  feed-line  vent  valve 

B.  CVC  Diffusion  pump 

1.  Close  2"  Te  me  seal  gate  valve 

2.  Turn  off  diffusion -pump  heater  switch 

3.  Let  water  and  Hypervac  25  run  until  heater  feels  cool 

to  the  touch 

4.  Turn  off  water 

•     •      5.     Turn  off  Hypervac  25 

C.  McLeod  pressure  system 

1.  Close  both  McLeod  gauge  shutoff  valves 

2.  Open  vent  line  valve 

3.  Turn  off  nitrogen  bottle 

4.  Open  needle  control  valve 
VII..    Periodic  Checks 

A.     Check  wooden  shims  under  the  test-tank  support  stand  twice 
weekly  and  knock  them  back  into  o  til.lons  .narked  on  the 
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floor  if  vibration  has  moved  them  out. 

B.  Check  the  oil  level  in  the  Hypervac  25  backing  pump  at 
least  every  other  day.     Add  Cenco  HyVac  oil  while  pump  is 
running. 

C.  Oil  level  in  the  Stokes  41 2H  pump,  should  be  halfway  up  the 
sight  glass  while  the  pump  is  running. 

D.  Temperature  of  diffusion-pump  cooling  water  taken  at  the 
outlet  should  be  110  -  120  degrees  Fahrenheit  for  maxi- 
mum pumping  rate.     Check  daily. 

E.  Stokes  pump  cooling  water  should  be  70  -  100  degrees 
Fahrenheit  at  the  outlet.     Check  during  each  period. of  ex- 
tended operation. 

F.  Drain  the  water  out  of  the  Stokes,  pump  and  refill  with  clean 
water  twice  weekly. 
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